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We study the dynamics of the collision between two fermions in Hubbard model with on-site 
interaction strength U. The exact solution shows that the scattering matrix for two-wavepacket 
collision is separable into two independent parts, operating on spatial and spin degrees of freedom, 
respectively. The S-matrix for spin configuration is equivalent to that of Heisenberg-type pulsed 
interaction with the strength depending on U and relative group velocity Vr- This can be applied to 
create distant EPR pair, through a collision process for two fermions with opposite spins in the case 
of \vr/U\ = 1, without the need for temporal control and measurement process. Multiple collision 
process for many particles is also discussed. 
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I. INTRODUCTION 

Pairing is the origin of many fascinating phenomena 
in nature, ranging from superconductivity to quantum 
teleportation. Owing to the rapid advance of exper¬ 
imental techniques, it has been possible both to pro¬ 
duce Cooper pairs of fermionic atoms and to observe 
the crossover between a Bose-Einstein condensate and 
a Bardeen-Cooper-Schrieffer superfluid [IHi. The dy¬ 
namic process of pair formation is of interest in both 
condensed matter physics and quantum information sci¬ 
ence. On one hand, the collective behavior of pairs gives 
rise to macroscopic properties in many-body physics. On 
the other hand, a single entangled pair is a promising 
quantum information resource for future quantum com¬ 
putation. 

In recent years, the controlled setting of ultracold 
fermionic atoms in optical lattices is regarded as a 
promising route to enabled quantitative experimental 
tests of theories of strongly interacting fermions M- 
In particular, fermions trapped in optical lattices can di¬ 
rectly simulate the physics of electrons in a crystalline 
solid, shedding light on novel physical phenomena in ma¬ 
terials with strong electron correlations [1, d, 0 • A major 
effort is devoted to simulate the Fermi-Hubbard model 
by using ultracold neutral atoms [lB - [i0 . This approach 
offers experimental access to a clean and highly flexible 
Fermi-Hubbard model with a unique set of observables 
pd| and therefore, motivate a large number of works on 
Mott insulator phase 0111 and transport properties 
[IE0, stimulating further theoretical and experimen¬ 
tal investigations on the dynamics of strongly interacting 
particles for the Fermi Hubbard model. 

In this paper, we study the dynamics of the collision 
between two fermions with various spin configurations. 
The particle-particle interaction is described by Hubbard 
model, which operates spatial and spin degrees of free¬ 
dom in a mixed manner. Based on the Bethe ansatz 
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solution, the time evolution of two fermonic wave pack¬ 
ets with identical size is analytically obtained. We find 
that the scattering matrix of the collision is separable 
into two independent parts, operating on spatial and spin 
degrees of freedom, respectively. The scattered two par¬ 
ticles exhibit dual features. The spatial part behaves 
as classical particles, swapping the momenta, while the 
spin part obeys the isotropic Heisenberg-type exchange 
coupling. The coupling strength depends on the Hub¬ 
bard on-site interaction and relative group velocity of two 
wavepackets. This finding can be applied to create distant 
EPR pair, through a collision process for two fermions 
with opposite spins without the need for temporal con¬ 
trol and measurement process. Multiple collision process 
for many particles is also discussed. 

The paper is organized as follow. In Sec. m we present 
the model Hamiltonian and analyze the symmetries. In 
Sec. mi we investigate some exact results obtained by 
Bethe ansatz concerning the two-particle problem. In 
Sec. mi we explore the dynamics of wavepacket colli¬ 
sion. Sec. |Vlis devoted to construct the scattering matrix 
for the collision process and the corresponding equivalent 
Hamiltonian. In Sec. ED we apply two-particle S-matrix 
to the case of multi-fermion collision. Finally, we give a 
summary and discussion in Sec. EhI 


II. MODEL HAMILTONIANS AND 
SYMMETRIES 

A one-dimensional Hubbard Hamiltonian on an A^-site 
ring reads 

N 

H = -K ^ (cl ^Cj+i^a + H.c.) + H (1) 

where c| g. is the creation operator of the fermion at the 
site i with spin a =t,i and U is the on-site interac¬ 
tion. The tunneling strength and the on-site interaction 
between bosons are denoted by k and U. For the sake 
of clarity and simplicity, we only consider odd-site sys- 
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tern with N = 2Nq + 1, and periodic boundary condition 

Ci,cr Ci_|_7V,cr' 

We analyze three symmetries of the Hamiltonian as fol¬ 
lowing, which is critical for achieving a two-particle solu¬ 
tion. The first is particle-number conservation [No-, H] = 
0, where iVo- = Si *^1 which ensures that one can 

solve the eigen problem in the invariant subspace with 
fixed TVo-j no matter U is real or complex. The second is 
the translational symmetry, [Ti,7J] = 0, where Ti is the 
shift operator defined as 


T, 






or T{ 
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t 

2 + 1 , 

-ik J 


( 2 ) 
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FIG. 1. (Color online) Schematic illustration of the collision 
process of two separated fermionic wavepackets with opposite 
spin orientations for three typical values of U. In all cases, 
the collisions result in momentum swap, but different spin 
conhgurations: (a) U — 0, swap; (b) U = oo, unchange; (c) 
U = ± |ui.|, maximal entanglement. 


This allows invariant subspace spanned by the eigenvec¬ 
tor of operator Ti. Based on this fact, one can reduce 
the two-particle problem to a single-particle probl 
The final is the SU(2) symmetry, [S^’^,H] = 0 
[5”^, =0, where the spin operators are defined as 


= 




t I ’ 


which satisfy the relation ] = 2S^. 






I vac), 




iKj 




(‘ 




± ^lA+r.. 


|vac), (r > 1), 


and 


and 

Ti (A)) = 14" (A)) , 

(11) 


Ti|<).± {K)) = {K)), 

(12) 

(4) 

ri+|+ (if))=e-^^^+|+ (A)), 

(13) 

while 

(5) 

+ 1^-(A)) = + !</>-(A)) = 0, 

(14) 


+ !<)>+ {K))=2\Spt +)), 

(15) 


++I+ (A)) =2+1+ (A)), 

(16) 

and 

ing 


5^ 1+(A)) = 5^ 1+ (A))=0, 

(17) 

(6) 

5-+I+ (A))=±+|+ (A)). 

(18) 


Then there are four invariant subspaces with {S,S^) = 
(0,0), (1,0), and (1,±1) involved. 


III. TWO-PARTICLE SOLUTIONS 


and 


T2 


1 giifr/2 y- 


piKjJ J 

S'.±t 


(r > 1), 


(7) 

where K is the momentum vector, indexing the subspace. 
These bases are eigenvectors of the operators fV^, Ti, 8“^ 
and . Straightforward algebra yields 


N,\c^-{K))=\cl>-{K)), (8) 

7 V,|<^±(iC))= |+(iC)), ( 9 ) 

iV^^|+(A)) = (l±l)E|^+(ir)), (10) 


Now we look at the two-particle solution in each in¬ 
variant subspace. We only focus on the solutions in sub¬ 
spaces (0, 0) and (1,0), since the one in subspace (1, ±1) 
can be obtained directly from that in subspace (1, 0) by 
operator S^. A two-particle state can be written as 

kl) = E (0 w {K )), (/^ (0) = /y, (-1) = o), 

r 

( 19 ) 

where the wave function (r) satisfies the Schrbdinger 
equations 


Qrfi,k (’’ + !) + Qr-ifK,k (’’-!) + 

[{-ir Qf Sr,No - Sk] f Ik ir)= 0 , ( 20 ) 
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FIG. 2. (Color online) Plots of the fidelity K^-tI {t))\ with 
the parameters Na = 20, = 62, ka = —ki, = tt/2, in the 
system Eq. © with = 81 and U = Vr- The red, blue and 
black lines represent the plots of fidelity |(^t| ^ (t))| in the 
condition of q = 0.13, 0.26, and 0.33, respectively. It shows 
that the fidelity is close 1, as a approaches to 0, which accords 
with the theoretical analysis in the text. 


and 


QrfK,k ir+l)+ Qr-lfK,k “ 1 ) + 

[Udr,0 + (- 1 )" Q^Sr,No ” (0 = 0 , ( 21 ) 

with the eigen energy sk in the invariant subspace in¬ 
dexed by K. Here factor = —2y/2KCOs{K/2) for 
r = 0 and —2kcos{K/2) for r 7 ^ 0, respectively. As 
pointed in Ref. [l 8 j | in previous works, the eigen problem 
of two-particle matrix can be reduced to the that of single 
particle. We see that the solution of ([^ is equivalent to 
that of the single-particle A^o -I- 1-site tight-binding chain 
system with nearest-neighbour (NN) hopping amplitude 
on-site potentials U and (—1)”^^ 2 k cos (Ai/2) at 0th 
and Noth sites, respectively. The solution of (l20l) corre¬ 
sponds to the same chain with infinite U. In this work, 
we only concern the scattering solution by the 0 th end. 
In this sense, can be obtained from the equivalent 
Hamiltonian 

00 

< = U |0) (0| + ^ (Qf K) (* + 1| + H.C.) . (22) 

i=l 

Based on the Bethe ansatz technique, the scattering 
solution can be expressed as 


where 

Aiy,/c = 4KCOs(Ar/2)sinfc, (25) 

A = 2tan~^ [ — ^ ) . (26) 

V ^K,k J 

And f, can be obtained from by taking U = 00 . 
We note that Rx^k {—U) = k = ^K-k, which reveals 
a dynamic symmetry of the Hubbard model with respect 
to the sign of U. 


IV. DYNAMICS OF WAVEPACKET COLLISION 

In this section, we investigate the dynamics of two- 
wavepackets collision based on the above formalism. We 
begin with our investigation from the time evolution of 
an initial state 


1$) = 1$,..) , (27) 

which represents two separable fermions a and 5, with 
spin a and cr', respectively. Here 

= |Vac), (28) 

with j = a, b and Na — Vf, ^ l/a, is a wavepacket with a 
width 2 -\/ln 2 /a, a central position N^ and a group veloc¬ 
ity u.,, = —2Ksmkj. We focus on the case (c, cr') = (t,4.)- 
The obtained result can be extended to other cases. In 
order to calculate the time evolution of state |$), two 
steps are necessary. At first, the projection of |$) on the 
basis sets {|(^+ (AT))} and {\(j)~ {K))} can be given by the 
decomposition 




Secondly, introducing the transformation 

Na=^iNa+Ni,),ra=Nb-Na, 

(30) 

kc — 2 ^ka ~t~ k}j^ , Qc — kf) ka-, 

(31) 

l=j + r, 

(32) 

and using the identities 

2 [(j - Naf + {l- 

= [{j +1)- {Na + N,)f + [(; - j) - (N, - 

A^a)]"(33) 

2 {kaj + kbl) 

= {ka + kb) {j + 1) + {kb - ka) {1 - j) , 

(34) 


fK,k ij) = + RK,ke^'"^, (23) have 


with eigen energy exik) = —4k cos (Ar/2) cos fc, k G 
[ 0 , 7 r]. Here the reflection amplitude 


RK,k 


i^K,k 't~U _ 
— U 


(24) 


= ^ V p-iK-2kcf/4a^ 

|^± (r„q,)), 


(35) 
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with 


where the S-matrix 


|^±) ^ ^^^g-a"(r-r,)V2gi9<=r/2 


(36) 


where ili 2 is the normalized factor. 

We note that the component of state |<i)) on each in¬ 
variant subspace indexed by K represents an incident 
wavepacket along the chain described by H^. This 
wavepacket has width 2-\/ln 2/a, central position Tc = 
Nb — Na and group velocity v = —4k cos {Kf2) sin (gc/2). 
Accordingly, the time evolution of state |<i)) can be 
derived by the dynamics of each sub wavepacket in 
each chain i^/g, which eventually can be obtained from 
Eq. (l23l) . According to the solution, the evolved state 
of {rcQc)) can be expressed approximately in the 

form e < 2 c /2 I'/’J (^cj ~9c)), which represents a 

reflected wavepacket. Here P (r/) is an overall phase, as 
a function of r', the position of the reflected wavepacket, 
being independent of U. In addition, it is easy to check 
out that, in the case with a 1, the initial state 
distribute mainly in the invariant subspace K = 2kc, 
where the wavepacket moves with the group velocity 
Uj. = —4 k cos (fcc) sin (gc/2) = Vb — Va- Then the state 
after collision has the approximate form 


1$ 


( 00 )) = -— 

|Vac)) 

I 

+^^e, 9 c /2 |Vac)) 

vSl ^ 

(^g-«=(;-JvO%*fc<.'4,^|Vac)), (37) 


which also represents two separable wavepackets at 
and respectively. Here is the normalized factor and 
an overall phase is neglected. We would like to point that 
the obtained conclusion is based on the fact that the 
shapes of two wavepackets |4)a_cr) and |4)f,_CT') are identi¬ 
cal. 


V. EQUIVALENT HEISENBERG COUPLING 


5 ^ g-*(e-7r)(^L-^R-i/4)^ ( 39 ) 

governs the spin part of the wave function. Here l.r 

denotes spin operator for the spins of particles at left or 
right, 9 = 2tan“^ [U/ (wr — ul)]^ where ul and ur rep¬ 
resent the group velocity of the left and right wavepacket, 
respectively. Together with the scattering matrix M for 
spatial degree of freedom 

M |L,p,ctl) |R,g,crR) = |L,g,crL) |R,p,crR), (40) 

we have a compact expression 

|f)=Af5|i), (41) 


to connect the initial and final states. In general the total 


scattering matrix has the form exp 


-iT Hdt 

J —00 


which 


is not separable into spatial and spin parts. Then Eq. 
(HU is only available for some specific initial states, e.g., 
spatially separable two-particle wavepackets with identi¬ 
cal size. This may lead to some interesting phenomena. 

It is interesting to note that the scattering matrix for 
spin is equivalent to the propagator 


S = T exp 


h (t) dt 


(42) 


for a pulsed Hensenberg model with Hamiltonian 

h(f) = J(f)(A-A-1/4), (43) 

with J J (t) dt = 9 — IT. Here T is time-ordered operator, 
which can be ignored since only the coupling strength 
J (t) is time dependent. This observation accords with 
the fact that, in the large positive U case, the Hubbard 
model scales on the t — J model [lanii , which also in¬ 
cludes the NN interaction term of isotropic Heisenberg 
type. 

This also indicates that the effect of collision on two 
spins is equivalent to that of time evolution operation un¬ 
der the Hamiltonian at an appropriate instant. 

In this sense, the collision process can be utilized to im¬ 
plement two-qubit gate. For two coupled-qubit system, 
the time evolution operator is simply given by 

W(t) = exp(-fl^L ■ ^Rt), (44) 


Now we try to express the two-fermion collision in a 
more compact form. We will employ an S-matrix to re¬ 
late the asymptotic spin states of the incoming to out- 
coming particles. We denote an incident single-particle 
wavepacket as the form of |A,p, ct), where A = L, R in¬ 
dicates the particle in the left and right of the collision 
zone, p the momentum, and cr =t, i the spin degree of 
freedom. In this context, we give the asymptotic expres¬ 
sion for the collision process as 

|L,p, ctl) |R, q, ctr) I— S |L, q, ctl) |R,P, ctr) , (38) 


yielding 

it) IDl IDr = e*‘/^cos i |t)L IDr - * sin ^ 11)^ |t)R), 

(45) 

where |cr =t,4-)L r denotes qubit state. We can see that 
at instants t = 7r/2 and tt, the evolved states become 

_iir/8 

^ (^/2) |t)L I'1 ')r = ^ (lt)L I'1 ')r “ * I'1')l lt)R) X46) 

WW|t)Lli)R = e-*"/"li)Llt)R, (47) 
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which indicates that U {tt/2) and U (tt) are entangling 
and swap operators, respectively. In practice, such pro¬ 
tocols require exact time control of the operation. 

Comparing operator U (t) and the S-matrix in Eq. 
(IMD, we find that two-qubit operations can be performed 
by the collision process, where U and relative group ve¬ 
locity Ur = ul — Ur, are connected to the evolution time 
by the relation 


t = 9 — TT = 2 cot ^ 

Then we can implement entangling and swap gates for 
two flying qubits via dynamic process. To demonstrate 
the result, we consider several typical cases with {7 = 0, 
00 , and ± |ur|, which correspond to the operations of 
swap, standby, and entanglement, respectively. The col¬ 
lision processes are illustrated schematically in Fig. [T] 
The advantage of such a scheme is that the temporal 
control is replaced by pre-engineered on-state interaction 
U. 

In order to check the above conclusion, numerical sim¬ 
ulation is performed. We define the initial and target 
states as 


14-(0)) = |L,p,t)|R, <?,!), (49) 

I^-t) = e*^(-zsin^ |L,g,t) \R,P,i) 

+ cos^ |L,g,4.) |R,p,t)), (50) 

where I'I't) possess the same relative position but the 
exchanged momentum compared with the state 14- (0)) 
as in Eq. (|4^ . On the other hand, we consider the 
evolved state 14- (t)) for the initial state being 14- (0)) 
driven by the Hamiltonian o, and caculate the fidelity 
I(4-7’! 4- (t))| in Fig. [2] It is shown that when the state 
14- (0)) evolves to the same position with |4'7’), the fidelity 
I (4-7-1 4- (t))| is almost to 0, which is in agreement with 
our previous theoretical analysis. 


VI. MULTIPLE COLLISION 

We apply our result to many-body system. Consider¬ 
ing the case that the initial state is consisted of many 
separable local particles with the same group velocity, 
termed as many-particle wavepacket train (MPWT), our 
result can be applicable if each collision time is exact 
known. In this paper, we only demonstrate this by a sim¬ 
ple example. We consider the collision of two MPWTs 
with particle numbers M and N {N ^ M). All the dis¬ 
tances between two adjacent particles in two trains are 
identical. The initial state is 


M N 

|LmjP: |Rni <Z: ’^n) j 

m—1 n—1 


OO 



t = oc 

Probability 



r M * 

i 

sin'(6»/2) 

- — 


* ♦ t 

i 

sin'(l9/2) cos' {6/2) 

11 * 


* t ♦ 

i 

sin'(0/2)cos'*(0/2) 



^ M t 

♦ 

cos^ {9 / 2) 



( * ^ 

♦ t 

1/8 



♦ * 

♦ t 

5/16 

* t 

( 

* ; 

* ♦ 

1/16 



« i 

* ♦ 

1/8 




* t 

5/16 



^ 11 


1/16 


(b) 

FIG. 3. (Color online) Schematic illnstration of the collision 
between the two MPWTs. (a) An incident single fermion 
comes from the left denoted as blue spin and collides with 
3-fermion train, which comes from the right denoted as red 
spins. It can be seen that the single fermion keep the original 
momentum, but it entangles with the 3-fermion train at the 
end of the collision. The amplitudes of the four states are 
listed. It is shown that the final state is direct product be¬ 
tween the states of single fermion and 3-fermion train when 
9 = TT, 8 = 0 with the corresponding parameter U = oo, 0, 
respectively, (b) The collision between the two MPWTs come 
from the opposite direction with particle number A = 2. And 
the probability for the superposition of states is listed with 
9 = 7r/2. 


where {Lm} and {R™} denote the sequences of particles, 
{am} and {t^} denote the spin configurations in each 
trains. According to the above analysis, after collisions 
the final state has the form of 

M N 

n n (52) 

m—1 n—1 

where the spin configurations {a'^} and {t(^} are deter¬ 
mined by the S-matrix, which is the time-ordered product 
of all two-particle S-matrices. During the collision pro¬ 
cess, the positions of particles in each train are always 
spaced by equal intervals. This makes it easier to deter¬ 
mine the times of each collisions. Then the final state 
can be written as 

M N M 

Si iLn, -1, Ai) |Rm:P; am) : (511) 

n—1 m—1 


where 


N 

s,=n 


Sl,N-n+l, 


n—1 


(51) 


(54) 
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and 




(55) 


where and k are corresponding Pauli matrices. Ap¬ 
plying the formula in Eq. (1551) to the case with M = 1, 
CTi =t, Tn =4-, n S [1, A^], we obtain 


N 


iLi,p,t) n 


■ i-7 • ^ (7 — 1) ^ 

—z > e sin - cos^-' ' — 
^22 
i=i 

X- ^ -i]^|L„,g,4.) |Ri,p,t) 


N 


-be*2^cos^ -]^|L„,(7,4.) |Ri,p,t) ■ 


(56) 


This conclusion is still true for the case with unequal¬ 
spaced {Rn}. For illustration, we sketch the case with 
M = 1, (Ti =t, N = 3, Tn =4-, n G [1,3] in Fig. |3l[a). 
One can see that the spin part of the final state is the 
superposition of the combinations of the four spins. 

Now we turn to investigate the entanglement between 
the single fermion and the MPWT with particle num¬ 
ber N. As is well known, the generation and control¬ 
lability of entanglement between distant quantum states 
have been at the heart of quantum information process¬ 
ing. Such as the applications in the emerging technolo¬ 
gies of quantum computing and quantum cryptography, 
as well as to realize quantum teleportation experimen¬ 
tally [12, Moreover, quantum entanglement is typ¬ 
ically fragile to practical noise. Every external manipu¬ 
lation inevitably induces noise in the system. This sug¬ 
gests a scheme based on the above mentioned collision 
process for generating the entanglement between a single 
fermion and the iV-fermion train without the need for 
the temporal control and measurement process. We note 
that although the incident single fermion keep the orig¬ 
inal momentum, it entangles with the 7V-fermion train 
after the collision, leading to a deterioration of its purity. 
To measure the entanglement between the single fermion 
and the A^-fermion train, we calculate the reduced den¬ 
sity matrix of the single spin 


Thus the purity of the single fermion can be expressed as 

P(oo) = Tr(p^)=2(A-l) +i, (59) 

where Tr(...) denotes the trace on the single fermion. For 
the case of A = 0, 1, we have P(oo) = 1, which requires 
0 — TT, and 0 = 0, obtained from interaction parameter 
U = oo, and 0, respectively. It indicates that the single 
fermion state and N-ieiimon train state are not entan¬ 
gled. In contrast, the purity P(oo) = 1/2 at A = 1/2 
when 

0 = 2cos“^ (2"5w) . (60) 

It corresponds to a completely mixed state of the out¬ 
going single fermion, or maximal entanglement between 
the single fermion state and A^-fermion train. Together 
with Eq. (IMl) . we have 


U = (UR 


ul) tan 



(61) 


which reduces to P ~ (ur — vi,)'/hi2N~^/‘^ for large N. 
This indicates that for large N, one needs to take a small 
U of order to generate the maximal entanglement 

between the single fermion state and A^-fermion train, or 
result in full decoherence of the single fermion. 

In the case of two-train collision, the calculation can 
still be performed in the similar way. However, it is hard 
to get analytical result for arbitrary system. Here, we 
sketch the case with M = 2, ai = <72 =t, A = 2, 
Ti = T 2 =4,, in Fig. mb). The probability on each spin 
configuration is listed as illustration. 


VII. SUMMARY 

In summary, we presented an analytical study for two- 
fermion dynamics in Hubbard model. We find that the 
scattering matrix of two-fermion collision is separable 
into two independent parts, operating on spatial and 
spin degrees of freedom, respectively, when two incident 
wavepackets have identical shapes. For two fermions with 
opposite spins, the collision process can create a distant 
EPR pair due to the resonance between the Hubbard in¬ 
teraction strength and the relative group velocity. The 
advantage of this scheme is without the need of temporal 
control and measurement process. Since it is now pos¬ 
sible to simulate the Hubbard model via cold fermionic 
atoms in optical lattice, these results can be realized ex¬ 
perimentally. 


Pnioo) = ^ ® ^ ^ , (57) 

where 

A = cos^'^^. (58) 
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